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The errors caused by qubit displacements from their prescribed locations in an 
ensemble of spin chains are estimated analytically and calculated numerically for 
a quantum computer based on phosphorus donors in silicon. We show that it is 
possible to polarize (initialize) the nuclear spins even with displaced qubits by using 
Controlled NOT gates between the electron and nuclear spins of the same phosphorus 
atom. However, a Controlled NOT gate between the displaced electron spins is 
implemented with large error because of the exponential dependence of exchange 
interaction constant on the distance between the qubits. If quantum computation 
is implemented on an ensemble of many spin chains, the errors can be small if the 
number of chains with displaced qubits is small. 
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I. INTRODUCTION 

A promising candidate for solid-state quantum computation is the phosphorus-doped 
silicon. This semiconductor material is the backbone of microelectronic technology. In 
order for this architecture to be useful for quantum information processing, the positions of 
the phosphorus dopants in silicon must be controlled by utilizing, for example, a scanning 



tunneling microscope (STM) , 0, y, 0]. The STM can be used to create many identical 
arrays of phosphorus atoms on the surface of silicon. A large number of such arrays is 
required for the detection of the qubit states. The long decoherence time of the nuclear spins 
of the phosphorus donors in silicon makes this quantum computer attractive for quantum 
information processing. 

Kane |6| proposed using nanoscale electronic gates to control the qubits. This technique 
has not yet been experimentally realized, so in this paper we consider a different architec- 
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ture. In our approach, the exchange interaction between qubits is constant, and selective 
interactions are reahzed through the use of a magnetic field gradient and both microwave and 
radio-frequency pulses. Measurement can be implemented using optical techniques similar 
to those used in Refs. [3, 0, 0] • 

Each phosphorus atom has a nuclear spin 1/2 and an electron spin 1/2. There is a 
hyperfine interaction between both spins. The information is stored in the states of the 
nuclear spins. The interaction between the nuclear spins of neighboring phosphorus atoms is 
mediated by electron spins coupled to each other by the exchange interaction. The advantage 
of an architecture based on controlled phosphorus impurities in silicon are the potential 
scalability and the possibility of using advanced silicon-based semiconductor technology in 
the quantum computer design. 

The typical size of an STM tip is larger than the lattice constant of silicon. Consequently, 
the phosphorus atoms in the lattice can be shifted from their prescribed locations by 1-4 
lattice sites. Errors are generated because the Larmor frequencies of the displaced spins and 
the exchange interaction constant are modified. In this paper, this error is estimated analyt- 
ically and calculated numerically for the protocol required to polarize the nuclear spins and 
for Controlled NOT gate between the electron spins of neighboring phosphorus atoms. If re- 
alized experimentally, each qubit chain in the first quantum silicon-based quantum computer 
will possibly contain the minimum number (two) of coupled qubits required to demonstrate 
basic principles of quantum computation in semiconductors. Hence, we specialize in this 
paper to the case of two phosphorus atoms in a chain. The unwanted displacements of the 
qubits seem to be inevitable in the quantum computer architecture. Therefore, our results 
are important for evaluating the possibility of constructing a working silicon-based quantum 
computer. 



II. HAMILTONIAN 

A schematic illustration of the system under consideration is given by Fig.^ If the qubits 
in each chain are placed at a separation of ~ 35 nm from each other, the nuclear-nuclear 
(~ 3.5 X 10~^ Hz), nuclear-electron (~ 0.65 Hz), and electron-electron (~ 1.2 kHz) dipole- 
dipole interactions are small compared to the electron-electron exchange interaction (~ 2 
MHz), so that one can neglect the dipole-dipole interactions. Since the relaxation time for 



3 



2^ = B^-Bf 

Hyperfine 
interactions 



t 



B 




X 

■4- 



Exchange 
interaction 



FIG. 1: A schematic illustration of two phosphorus atoms placed in a permanent magnetic field. 
The electron spins Si and «S'2 (large arrows) of the neighboring atoms interact with each other via 
the exchange interaction and the nuclear spins Ii and I2 (small arrows) interact with the electron 
spins through the hyperfine interactions. 

the electron spins at temperatures of 1-7 K is relatively short (0.6-60 ms at 7 K Q), the 
quantum information must be stored in the states of the nuclear spins. Because the nuclear 
spins do not interact directly, electron spins can be used to mediate the nuclear-nuclear 
interactions. In this setup, the electron spins must be coherent only during the relatively 
short time of implementation of the quantum logic gate, such as a Controlled-NOT gate 
on a particular pair of qubits. The coupled electron spins of the neighboring phosphorus 
atoms serve as the static "quantum bus" which transfers the interaction between nearest- 
neighbor and distant nuclear spins in a chain with more than two phosphorus atoms. In a 
similar technique "flying" qubits represented by photons has been proposed to mediate the 
qubit-qubit interaction between trapped ions |ll[ [l^, Q, Ili > 15. Ila Irn]. 
The Hamiltonian is 

H = H'^ + V{t), (1) 

where is the unperturbed Hamiltonian of the system shown in Fig.Cland V(t) is the time- 
dependent field describing rectangular radio-frequency pulses used to implement quantum 
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logic gates. The Hamiltonian H'^ is 

= ^,BISI + ^eB^Sl - ^nBfn - inBfn + A ( SJ\ + S2I2 ] + (2) 



where Ik and 5*^ are the spin operators of, respectively, the fcth nuclear and electron spins, 
k = 1,2; /| and Sf. are the projections of these operators on the z axis; 7e and 7„ are, 
respectively, the magnitudes of the electron and nuclear gyromagnetic ratios; B^ is the 
permanent magnetic field at the location of the kth spin; A and J are, respectively, the 
hyperfine and exchange interaction constants. 

We utilized perturbation theory to calculate the eigenvalues Ei, i = 0, 1, . . . , 15 of the 
Hamiltonian which are given in Appendix In spite of the fact that our system 

allows exact analytical calculation of Ei, we prefer to use a perturbation approach because 
this technique allows us to separate the most important contributions E^^^ to Ei from the 

less important and to use relatively simple expressions for the eigenvalues to analyze the 

(2) 

dynamical properties of this system. The small corrections E^ are useful for calculating 
the parameters of the pulses in numerical simulations. The eigenvalues E^, i = 1, 2 . . . , 14 
are calculated with an accuracy of ^^(A/2)/(27r) ^ 23 Hz, where 

A 

e = r ~ 6 X 10"^ 

Here we assume b = (5° + B^)/2 = 3.3 T, so that 7e6/(27r) = 92.5 GHz, 7e/(27r) = 28.025 
GHz/T, A/{2'k) = 117.53 MHz Q- 
The basis states are 

\n2e2e1n1), (3) 

where the electron spin and nuclear spin rii of the zth phosphorus atom, i = 1,2, can 
assume the values and 1; the state |0) corresponds to the orientation of the spin along 
the direction of the permanent external magnetic field and the state |1) corresponds to the 
opposite direction. The state |0000) (which is not the ground state) has energy Eq, the state 
|0001) has energy Ei, etc. In the spin chain Q there are interactions only between the 
neighboring spins, so this kind of spin ordering is convenient for the analysis of conditional 
quantum logic gates. 

The basis states are not the eigenstates of the Hamiltonian Hq because of the (off- 
diagonal) terms AI^S%, AljiSl, k = 1,2, and JS'fS'f, J^fS*!. However, the eigenfunctions 
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approximately coincide with the basis states Q if the conditions 



27e(5| - Bl) ' |27e(5| - Bl) - A 

are satisfied ^18]. The unwanted changes (errors) in the wave function due to the infiuence 
of the off-diagonal terms are of the order of e or e' or ^. 



III. IMPLEMENTATION OF QUANTUM LOGIC GATES 



The time-dependent magnetic field has the following components: 



B^(t) = B\cos{ut + ip),- sm{ut + <^), 0), 



(5) 



where B^, u, and ip are, respectively, the amplitude, frequency and phase of the pulse, 
and t is time. The frequency u can assume both positive and negative values The 
time-dependent term in the Hamiltonian has the form 



{St 



\ [it + It 



g-i(i/t+v.)_l_ 



(6) 



where n° = -f^B^ and = -fnB\ 

We will consider the following three basic kinds of quantum gates: a Controlled NOT 
gate with the control nuclear spin and the target electron spin of the same phosphorus atom, 
a Controlled NOT gate with the control electron spin and the target nuclear spin, and a 
Controlled NOT gate between the neighboring electron spins. We now briefiy describe the 
procedure for calculating the optimal parameters of the pulses required to implement these 
gates. (For a detailed analysis see Ref. ^^.) Assume that the direction of the kth spin 
in the state \p) is along the direction of the permanent magnetic field -B°, k = 1,2 (i.e. 
\p) = I ■ ■ ■ Ofc ■ ■ ■)) and the state |g) = \ ■ ■ - Ik - ■ ■) is related to the state \p) by a fiip of the 
A;th spin. In order to fiip the kth spin in the state \p) or in the state the frequency u 
and the time-duration r of the rectangular pulse must satisfy the conditions 



Eg — Ep, 



^ = -, 

TT 



(7) 



where fl is the Rabi frequency of the pulse. 
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FIG. 2: A schematic illustration of the energy levels ivolved in implementation of a Controlled-NOT 
gate. The frequency v is the frequency of the pulse. 

Assume that we wish to suppress the transition \p') ^ |g'), where the fcth spin in the state 
is along the direction of the permanent magnetic field B\ and the state is related 
to the state by a flip of the /cth spin; \p) ^ and \q) ^ \q'). A schematic illustration 
of the involved energy levels and transitions is shown in Fig. |21 The Rabi frequency of the 



pulse must satisfy the 27r_ft'-condition 
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2l| 



|A| 



(8) 



where i^' = 1, 2, . 
nuclear spin, and 



is an integer number, Vt = Vl^. for an electron spin and Vt = Qn for a 



A = Eg, - Ep, - u. 



(9) 



The probability of transition generated by the pulse with the Rabi frequency fl and detuning 



A is ^ 



R 



ttA 



A^"" 2?^ 



(10) 



where A = v/TFTA^. We have i? = 1 for A = and i? = if satisfies Eq. These 
two conditions allow one to implement the Controlled NOT gate as indicated in Fig. El 



IV. J AS A FUNCTION OF THE DISTANCE BETWEEN THE QUBITS 



A simple approximation of t 
by the Herring-Flicker formula 



le exchange interaction adequate for our purposes is given 

2 / ^ N 5/2 



J = 1.642- 



a y 



2a\ 
exp ( 



(11) 
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where e is the electron charge, k = 11.9 is the dielectric constant of silicon, and a is the 
distance between the qubits. The effective Bohr radius is 

aB = K—al ^ 33.14 A, (12) 

where 0.5292 A is the Bohr radius, M is the electron mass, M* = 0.19M is the effective 
electron mass. In Eq. ()11|) . a = Na^ where — 1 is the number of interstitial silicon atoms 
between the two phosphorus atoms serving as qubits and Oq = 7.68 A. We assume that the 
qubits are placed on the (100) surface of silicon and the direction of the permanent magnetic 
field is perpendicular to this surface. In Table I we show the values of J/ (27r) for different 
calculated using Eq. (fTTjl . The controlled separation between the qubits can be varied by 
changing N. In this paper we assume that the predefined spacing (desired spacing without 
unwanted displacements) is Nq = 47, so that the predefined exchange interaction constant 
is Jo/(2vr) = 1.97 MHz. 



N 
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41 
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48 


49 


50 


51 


a (nm) 


30.72 


31.49 


32.26 


33.02 


33.79 


34.56 


35.33 


36.10 


36.86 


37.63 


38.40 


39.17 


J/(27r) (MHz) 


33.75 


22.58 


15.09 


10.07 


6.71 


4.465 


2.97 


1.97 


1.306 


0.865 


0.573 


0.37855 



TABLE I: The distance a between the neighboring phosphorus atoms and the exchange interaction 
constant J/{2ir) for different N, where — 1 is the number of interstitial silicon atoms. 



V. CONTROLLED NOT GATES BETWEEN THE ELECTRON AND NUCLEAR 

SPINS 

We now derive the pulse parameters of four pulses required to initialize the nuclear spins. 
Each pulse implements one logical gate. The gates are jl8| 

(a) CNni,el, (b) CNel,„l, (c) CN„2,e2, (d) CNe2,n2, (13) 

where CNj ,,- denotes a Controlled NOT gate, i is the number of the control qubit and j is 
number of the target qubit. The protocol (jl3p includes also a delay time (of the order of 0.1 
ms) between pulses (b) and (c) and after the pulse (d). During this time the electron spins 
relax to their ground states. We assume that the relaxing electron spins do not influence 
the corresponding nuclear spins jlS^ . 
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In order to calculate the energies and resonant frequencies we will use the following 
parameters: 

^ = 0.985 MHz, 4^ = 58.765 MHz, ^ = 92.48 GHz, ^ = 56.93 MHz, 
2tt ztt ztt ztt 

^ = 65.76 MHz, ^ = 40.48 kHz, (14) 
2tt 2tt 

where 6B = {B2 - Bi)/2. Here we assume h = 3.3 T, 7„/(27r) = 17.25144 MHz/T, the 
magnetic field gradient is 1.3 x 10^ T/m, and from Table I we use a = 36.1 nm. 

We now calculate the parameters of the pulse implementing gate (a). The resonant 
transition is |1101) ^ |1111). The resonant frequency for this transition is 



_ Ei5 - ^{3^ _ 1 
27T ^ 2tI 271 



A J 



-leh + - + - + ^{le5Bf + .P/A 



-92.35 GHz. (15) 



The near-resonant transition to be suppressed is |1100) ^ |1110). The detuning for this 
transition is 

A . E 



^1 -£^14 ~ -£^12 ~ 



27r 2ti 



1 

2^ 



(7e5i? - A/2Y + JV4 - ^{ie5BY + J74 - ^ 



-117.47 MHz. (16) 



In order to suppress the near-resonant transition with the detuning Aj the Rabi frequency 
of the pulse must satisfy the 2TiK condition 

I A,- 1 



(17) 



277 2^^AK'^ - 1' 

where K is an integer number and i is the pulse number. For the first pulse with Ai given 
by Eq. (fTBj) and ioi K = K^, = 1, we have 

— = = 67.82 MHz. (18) 

27r 27rv3 

There are three sources of errors in our model: (i) nonresonant excitations, (ii) off- 
diagonal components of the exchange and hyperfine interactions, and (iii) qubit displace- 
ments. The error in the probability amplitude due to the nonresonant transitions is propor- 
tional to the ratio j^l 

where 5uj is the difference between the frequency corresponding to the nonresonant transition 
and the frequency of the pulse. 
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FIG. 3: A schematic illustration of an ideal qubit chain (upper part) and a qubit chain with a 
displaced qubit (lower part). 



The order of magnitude of the error generated by a single pulse acting on an electron 
spin can be estimated by adding probability errors generated due to all three mechanisms 



1 - m 



P = 1- R + e^ + I?, 



(20) 



where e, i?, and n are given by Eqs. (jH), (fTUj). and (fT^. 

We now analytically estimate the error generated by the first pulse due to the displace- 
ment of the first phosphorus atom using Eq. (llUj). This displacement introduces an unwanted 
change in the Larmor frequency of the first qubit and the correction dJ to the exchange in- 
teraction constant J. Let qubit 1 be displaced in the negative x direction by one lattice site 
as shown in Fig. El From Table I we find 

dJ 



27r 



1.306 MHz - 1.97 MHz = -0.664 MHz. 



(21) 



The permanent magnetic field at the location of the first qubits is B\ — dB, where 



From Eq. (fT3)) we have 



- B^ 

dB = i ^ 9.985 X 10"^ T. 



^ -7,5° + Y + 2 ' 



(22) 
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FIG. 4: Probability error Pi generated by gate (a) in Eq. as a function of displacement m of 
the first phosphorus atom for four values of Kg. 



so that the unwanted detuning is 

A[{m = -1) 



27r 



2.466 MHz, 



(23) 



where 'yJB/{27T) = 2.798 MHz and we used Eqs. (ED) and If we put instead of A 
and Qi instead of Q in Eq. ()10|) we find the probability error 1 — R due to the unwanted 
qubit displacement. 

For displacement in the opposite direction with m = 1 



1 

2^ 



5 J'' 



-1.8 MHz, 



(24) 



27r 27r V ' 2 ) 

where 5J' /{2-k) = 1 MHz (see Table I). Since \A[{m = 1)| < \A[{m = -1)|, the error 
Pi for the displacement in the negative direction along the x axis is larger than the error 
for displacement in the positive direction because 1 — i? in Eq. (p]]) increases with \A'i\ 
increasing. 

In Fig. |3]we plot the error Pi for the transition |1101) — ^ generated by the first 

pulse for the chain with displaced qubit 1 as a function of the qubit displacement m along 
the X axis. This plot is asymmetric with respect to m = which follows from the analysis 
presented in the previous paragraph. The data are obtained by utilizing the eigenfunctions 
of full Hamiltonian in the rotating frame. (For a description of the numerical solution 
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FIG. 5: Probability error Pi generated by by gate (a) in Eq. ()13() as a function of displacement m 
of the first phosphorus atom along the x axis for K^. = 2. The estimate is calculated using Eq. (j^U)) . 

see Ref. The parameters of the pulses in our numerical simulations are calculated 

analytically using the eigenvalues from the Appendix. In the figure, m is the number of 
lattice sites for a displaced qubit, so that the displacement is mao, m = —4, —3, . . . , 4. Our 
simulations (not presented here) show that this and other pulses generate the maximum 
error in situation when the "active" qubit to be flipped is displaced and the error is smaller 
when the neighboring qubit is displaced. 

Since \A[\ ^ Qi [see Eqs. (fTK|l . (j2SI), and (j21I)] the error Pi is relatively small small. 
(Here and below we assume that the error P is small if P ^ 1.) The smallness of the error 
generated by the first pulse is due to (a) the neighboring phosphorus atoms being located 
relatively far from each other; (b) the electron Rabi frequency being relatively large, and (c) 
the exchange interaction constant A, responsible for the Controlled NOT gate, not being 
modified by the displacements of the phosphorus atoms. 

From Fig. m one can see that the error increases as \m\ increases. As follows from this 
figure, the probability error becomes more sensitive to the qubit displacement when Ke 
becomes relatively large and the Rabi frequency becomes relatively small. In Fig. El the 
probability error estimate ()20j) is compared with the results of the numerical solution. One 
can see that our estimate is satisfactory and supports the approach taken. The asymmetry 
of the plot in Fig. with respect to m = is of the same origin as that in Fig. 0] discussed 
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before. 

We now estimate the error associated with a flip of the flrst nuclear spin in the situation 
when the flrst phosphorus atom is shifted from its ideal position to m = —1. We flrst 
calculate the pulse parameters. The resonant transition is ^ |1110). The frequency 

of this transition is 

^ - ^11^ = ^ Lb - ^JB + ^)= 115.655 MHz. (25) 

The near-resonant transition to be suppressed is |1100) ^ |1101). One can show that the 
detuning for this transition is A2 = Ai, where Ai is given by Eq. ()16p. The 2'kK condition 
for the second pulse reads 

0„ IA„I 

(26) 



'2 



271^ AKl - 1 

where is an integer. In order to suppress the flip of nuclear spin 2, the value of ^2 must 
be much smaller than the frequency difference associated with flip of the second nuclear spin. 
The unwanted off-resonant transition is |0110) ^ |1110). The detuning for this transition 
is 5uj = 2'jnSB. From Eq. (fT^ and for /i <^ 1, we obtain the condition 

> 363, (27) 

or ^2/(271) < 162 kHz. 

We now calculate the error associated with displacement of the flrst phosphorus atom 
for m = —1. The corrections to the magnetic fleld is given by Eq. (j22j) . The correction 
to 7„5-B/(27r) is 0.86 kHz and the correction to 7„6/(27r) is -0.86 kHz. From Eq. (j25p the 
unwanted detuning from the resonance is 

A'o 

= -1.72 kHz. 

2tt 

From the condition Q2 ^ I '^2!) '^^ S^^d Q2 ^ 1-3 kHz, or 

< 34, 148. (28) 

From Eqs. ()27|) and ()28|). one can see that it is possible to make both the errors due to the 
nonresonant excitations and the errors due to the unwanted displacements small. 

In Fig.iniwe plot the error as a function of dimensionless displacement m for four different 
values of Kn- This plot has the same features as Fig. IHa). The probabihty error P2 is 
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FIG. 6: The probability error P2 generated by the second pulse on the first nuclear spin as a 
function of displacement m of the first phosphorus atom along the x axis for five values of K^- 
The data are obtained by numerical solution of the full Hamiltonian . 

independent of the direction of the displacement because P2 depends only on the magnitude 
of the deviation of magnetic field from the optimal value in the location of the qubit. The 
error due to the deviation of J from the optimal value Jo does not contribute to P2 because 
the first nuclear spin does not directly interact with the second electron spin and second 
nuclear spin. Since the nuclear spins do not interact with each other and electron and nuclear 
spins of different phosphorus atoms do not interact with each other, the displacement of the 
second phosphorus atom practically does not affect the quantum logic operations on the 
nuclear spin of the first atom. This was also confirmed by our calculations (not presented 
here). This result follows also from Eq. where z/2 = 7n-Bi + A/2 is independent of J 
and B^. 



VI. ERROR IN THE INITIALIZATION ALGORITHM 

We have analyzed the errors generated by the first and second pulses. The errors gener- 
ated by the third and fourth pulses [pulses (c) and (d) in Eq. are of the same order 
of magnitude. As follows from these results, one can make the error relatively small and 
polarize (initialize) the nuclear spins even in the situation when the phosphorus atoms are 
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Kn 

FIG. 7: The probability of error P = 1 — |Ci5(T)p generated by the initiahzation protocol as a 
function of Kn- Each phosphorus atom in an ensemble of 10,000 spin pairs had the probability rm, 
m = —4, —3, . . . , 3, 4 to be displaced in a random direction along the x axis. For a perfect sample 
without displacements = (solid line), Kg = 1. The data are obtained by numerical solution 
of the full Hamiltonian 

displaced from their prescribed locations in the chains. To illustrate this feature we modeled 
the initialization protocol on an ensemble of 10,000 noninteracting spin chains. Each chain 
consisted of two phosphorus atoms. The initial state for each chain was the superposition 

C6(0)|0110) + C7(0)|0111) + Ci4(0)|1110) + Ci5(0)|llll) 

with randomly chosen normalized complex coefficients Cj(t), i = 6,7, 14, 15. The probability 
of error, P = 1 — \Ci^{T)\'^, where T is the time of implementation of the initialization 
algorithm, as a function of Kn is presented in Fig. [7| We calculated three situations. In the 
first situation each phosphorus atom in the ensemble had the probability = (1/2) 1™'"''^, 
m = —4, —3, . . . , 3, 4 to be displaced in a random direction along the x axis (filled circles 
in Fig. [7j). In the second situation the probability of displacement is = (1/2)!™'+^ which 
corresponds to a sample of higher quality. For comparison, for the third situation we plot P 
for a perfect sample without displacements (r^ = 0). The data in Fig. [7| are averaged over 
40 realizations with different initial conditions and different sets of random displacements. 
The size of the error bars is of the order of the size of the symbols. The Rabi frequencies 
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corresponding to different values of Kn can be obtained from Eq. ()26|) witli |A2|/(27r) ^ 
|A4|/(27r) 117.5 MHz. If Kn assumes the values in the interval from 700 to 10,000, the 
magnitudes of ^2 and vary in the interval from 84 kHz to 5.9 kHz. 

As follows from Fig. 6, when Kn is relatively small (and the Rabi frequency is relatively 
large) the error is mostly determined by the off-resonant excitations and only slightly depends 
on the quality of the sample. For relatively large values of Kn (and small Rabi frequencies) 
the error is smaller for better samples with smaller number of displacements. 



VII. CONTROLLED NOT GATE BETWEEN ELECTRON SPINS 

As follows from the results presented above, one can make the error in the Controlled 
NOT gates between the nuclear and electron spins of the same phosphorus to be relatively 
small because the coupling constant A responsible for the Controlled NOT gates is not 
modified by the displacements. Now we will consider the error generated by the Controlled 
NOT gate between the electron spins of the neighboring phosphorus atoms. One can expect 
this error to be relatively large because a displacement modifies the coupling constant J 
which is the key parameter for this gate. 

The detuning is |Ae| ^ J , so that the Rabi frequency is equal to 

= , ^ , (29) 

^|AKf^l V^4irf - 1 

where K'^ = 1,2,... is an integer number. From Eq. (fTT|) we can find the modification 
5J = J — Jq from the "ideal" value Jo due to the displacement aom = 7.68m A, 



^ I „2aom/aB _ i 

2N^ 



For = 47 and = 33.14 A, we obtain 

6 J/ Jo = 0.4103m + 0.1082m2 + 0.0166m=^ + 0.0019m^ + . . . . 

For the minimum possible displacement with m = —1, we have SJ/Jq ~ —0.32 which can 
be verified from data in Table I. 

The unwanted modification of the Rabi frequency ()29j) caused by the displacement is 
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The error is small if two conditions are satisfied: (i) \SQe\ ^ and (ii) modulus of the 
unwanted detuning |Ag| is small, |Ag| <^ Q^.- One can show that the unwanted detuning 
from the resonance caused by the displacement is |Ag| ^ 6J/2, so that both conditions (i) 
and (ii) cannot be satisfied by any choice of parameters of the system and the error is large. 
We checked this result numerically and determined that the error is close to unity for 
|m| ^ and \Pe\ ~ 10"^ to lO'^ for m = 0. 

We now estimate the error if one qubit is displaced in a spin chain with more than two 
phosphorus atoms. The presented above analysis is valid for the edge qubits of this chain. 
Consider, for example, the state | ■ ■ ■ Ik+ilk^k-i ■ ■ ■) formed by coupled electron spins. (Here 
we do not indicate the states of the nuclear spins.) If the the kth atom is not displaced, 
then the frequency resonant for the transition of the kth qubit is 

27r 27r V 2 2 / 27r ' ^ ^ 

where Jq/{2tx) = 1.97 MHz. If the kth qubit is displaced in any direction along the chain by 
one lattice cite then instead of the term Jq/{2ti) in the right-hand side of Eq. (|!?T|) we have 
from Table I 



1 

2^ 



J(m = 1) Jijn = — 1) 



2 2 

so that the unwanted detuning from the resonance is 



2.14 MHz, 



A'g = O.O86J0. (32) 

From Eq. ()29|) for Kg = 1, we obtain 

n^ = ^^ 0.58Jo (33) 
v3 

From Eqs. ()32j] and ()33j] one can see that the unwanted detuning is relatively small | Ag| '^flg 
and condition (ii) is satisfied. One can show that this condition is not satisfied for the 
transitions | ■ ■ -Ife+ilfcOfc-i ■ ■ ■) |lfc+iOjfcOfc-i). The deviation in Rabi frequency [condition 
(i)] is always relatively large. For example, if the target kth qubit is displaced in the negative 
direction (m=-l) and the {k — l)th qubit is the control qubit, then 6fle is given by Eq. (IHU)) 
and condition \6Qe\ ^ is not satisfied, so that the Controlled NOT gate between the 
displaced electron spins generate large error. 

Presented here analysis demonstrate that the Controlled NOT gate between the electron 
spins cannot be performed on displaced qubits. If quantum computation is implemented 
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on an ensemble of many spin chains, the errors can be small if the number of chains with 
displaced qubits is small. 

In conclusion, we note that the problem of unwanted displacements appears only in 
quantum computation performed on an ensemble of spin chains. If one works with one chain 
only, one can measure the positions of the phosphorus atoms using a scanning tunneling 
microscope and optimize the parameters of the pulses appropriately. 
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APPENDIX 



The eigenvalues of the system of two phosphorus atoms in a permanent magnetic field 



are 



^ / ^. ^ J 

Eq = (7e -7n)0+ - + -, 



4 4 



J 



.(2) 



^2 



E?^ - Ef^ 



J 



J2 



^(2) 



.(2) 



4 
J 



e: 



(2) 



= -7n&-T-\/(7e5S)2 + :-, El 



J 



J2 



.(2) 



3 -0, 



p(0) p(0) 



(0) 



5 4 



.Er 



E, 



(0) ' p(0) 



Er - E, 



(0) 



(0) 

'6 



-(7e + 7n)6- 2 + 4' ^6 =T 



Ef' - Ef' 



+ 



p(o) p(0) 
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E^'^ ^-jeb-7n5B + ^, Ef^ 



p(o) p(0) 



E^^^^^b + ^JB + i Ef^^-E?\ 



E, 



(0) 



(7e + 7n)&- 2 + 4' ^9 =T 



1 



Ef - Et 



+ 



1 



p(o) _ p(0) 



^10 - 4 

J 



p(0) p(0) 



+ 



p(o) _ p(0) 



J2 



.(2) 



.(2) 



e[? = 7„5S - 



4 
J 



J 



J2 



(0) 
14 





J2 

4 ' 


P(2) - 


= 0, 


p(2) 
-C'13 








4 


rp(o) 

-^13 


p(0)l 


p(2) 


p(2) 
- -J^13 , 





-7e^ + 7n55 + ^, 
^15 = (-7e + 7n)& + + ^' 

where b = {B^ + B^)/2 and 55 = (5^ - Bf)/2. In the text we assume Ei = E^ + £;f \ 
The eigenvalues £^io and £^12 are presented for the case A/2 < je^B. For the opposite case, 
A/2 > 7e5S, one must exchange the eigenvalues £'10 <-> -E'i2- 
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